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Control design plays a major part in electrostatic linear comb actuator (ELCA) 
systems. Most existing control methods of ECLA mainly focus on traditional 
PID method. However, proportional—integral—derivative (PID) controller is 
only appropriate for the reference to be known as constant. In this article, the 
problem of establishing the dynamic model and designing the controller is 


investigated for ECLA working in the time-varying reference. The tracking 
effectiveness is constructed for ECLA systems based on the consideration of 
roots of the equivalent differential equation in tracking error model. 
Moreover, the convergence velocity is also tuned by changing the parameters 
in the presented control scheme. Simulation results guarantee the performance 
of the nonlinear control scheme as well as dynamic model to be established. 
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1. INTRODUCTION 

Micro electro-mechanical system (MEMS) is a remarkable device and popularly considered in several 
directions such as robots, vehicles, aerospace, industrial electronics, medicine, etc. Micro actuators can be known 
as one of the MEMS devices as well as describe an efficient task creating a force to move the micro elements and 
micro systems. Electrostatic linear comb actuator (ELCA) is a remarkable instrument structure using the working 
description via electrostatic force for establishing linear motion. The ELCA systems have many advantages 
including high effectiveness, low power requirement, and simple structure. They are able to explain the fact of 
easy movement in MEMS elements, such as resonator systems [1], gyroscopes [2] and micro systems of the 
conveyer or transportation systems [3], micro motors or gripper [4], [5]. Several work of ELCA have been 
recently developed with many appropriate results including new structures, working effectiveness, 
manufacturing improvement, smart materials. Particularly, the fabrication effectiveness and the ELCA’s 
movement stability guarantees the operating effectiveness of the elements. The characteristic of suspension 
beams with a crab-leg, a fold [6], [7] was developed to decrease the stiffness and change the deviation. A 
different direction on improving deviation was addressed [8], by decreasing the movement of out-of-plane 
from that generates the in-plane deviation. The work in [9] presented an actuator property after inserting more 
a module of a sensor as well as a sub-actuator for stability of control systems. The effect of side etching on 
deviation of ELCA with disability rent increase exemption DRIE method was pointed out in [10]. By 
computing the ELCA deviation in design objective [3], [4], [10]-[12], it is necessary to utilize the theoretical 
method to be developed by a static equilibrium equation and decreases a design time. 
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The fact is that these elements manipulate in statement of motion by using AC voltage such as 
traditional waves. Hence, under the effect of force and damping follows an efficient property in estimating the 
deviation of ELCA, which was considered in [13], [14] with the EAs to be developed for moving resonators, 
gyroscopes. Moreover, the motion dynamic model including an efficient mass and factor of conversion of air 
damping was not considered in these references. To achieve better behavior, the consideration of using control 
scheme is becoming a promising method [8]. However, the controller of ECLA have been considered with the 
traditional PID method. It should be noted that conventional PID controller does not satisfy the tracking 
effectiveness under the time-varying reference. Moreover, nonlinear controller and optimal control scheme 
have been usually proposed for different systems, such as robotic systems, power electronics [15]-[25]. 
Motivated by the above-mentioned discussions, this article studies the dynamic model of ECLA and the 
corresponding controller. To improve the tracking effectiveness under time-varying reference, the controller is 
carried out based on the consideration of roots of differential equations. 


2. PROBLEM STATEMENT 

Figure 1 describes the structure of ELCA including a moving element © and fixed element (@ and ©). 
The moving element with moving comb fingers is pressured by four EB systems driving in y-direction based on 
the electrostatic force (EF), one terminal of each beam is kept on ©. The configuration of movable and fixed 
comb fingers can be seen with the same shape. The work in [4] shows the ECLA geometrical parameters in Figure 
1, where b is a beam width, A is the device layer thickness, and gz, is the gap. 
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Figure 1. Structure of ELCA 
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Under the influence of a voltage V between © and © of the ELCA, the total force in y-direction can 
be computed as (1). 


F, = Nhe EG y2 (1) 
Jo 


Where £o = 8.854 x 10°°pF/um and s= 1; V is the moving voltage. The stiffness in y-direction with four 
suspension beams system is described as (2) [4], [10]. 


_ 4Ehb? 
aaa 


(2) 


Where E = 169GPa. 
Due to the application of a DC voltage for ELCA, the deviation yo of ELCA (i.e. displacement of the 
shuttle) is directly computed from the equation of force balance: 


h 4Ehb? 
F, — Fa = Oor: =v? = Mo 
It implies that: 
— nééq (L 3 2 
n= y o 
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In (3) shows the quadratic relation between the deviation of the ELCA and the DC moving voltage. The 
objective of this article is to establish the dynamic model of ELCA for designing the controller with the 
reference of output y(t) depending on time. 


3. DYNAMIC MODEL OF ELCA AND CONTROL DESIGN 

Under the AC moving voltage, it implies the movement of ELCA’s shuttle to be similar with the 
movement to be described in Figure 2. Where K, C, M are the coefficients of conversion stiffness, 
conversion damping and conversion effective mass, respectively; F,(t) is an electrostatic moving force. 


Fe(t) 


Figure 2. 1-DOF dynamic model of the ELCA 


The dynamic model in y-direction (Figure 1) is represented as (4). 
My + Cy + Ky = R(t) (4) 


For the requirement of calculating exactly the deviation of the ELCA using dynamic model, the factors K, M, 
C are necessary to be accurately identified. Moreover, the factor K in y-direction is given in (2) and the 
appropriate mass of movable elements, including the elastic beams, is transformed into shuttle’s motion with 
the principle of the kinetic energy of mass M being equal to all movable parts’s kinetic energy. It follows that: 


í 2 
M =M, +f" 22) dx (5) 


“EL 20 Ga 
moreover, the related equation between y(x) and y(L) is given through the relation of deviations: 


3Lx? -2x3 


2- 3 . . 
yw = yh) > 9) = FD (6) 
according to (5) and (6), it leads: 
13M 52bL 
M=M,+4 se = phh + =) (7) 


where p = 2330 kg/m? is the silicon mass density; A; is the total top side area of the shuttle. Based on [15], the 
coefficient C in (4) can be calculated: 


C=C, +C +C 


where 
c = u (+s) (8) 
C, = ys (9) 
Cs x Sut (10) 


combination of in (8), (9) and (10), we have: 


L (Abt 4 2Ase | AttAss | 16h 
a e Eae re) (11) 


considering the driving voltage with the square wave (Figure 3). 
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Figure 3. Square wave of driving voltage 


The force function F. (t) is represented as: 


NnhEEo 172 
Vé O<t<T/2 
RO= o" OSES a2) 
0 T/2<t<T 


Yı = eP! (A, cos( wt) + Az sin(wt)) + Ë (13) 
Here: B = £, w = 4J w2 — B; 0, = IE With the initial condition: y(0) = Oand y(0) = 0, it follows: A, = 
= fe, A, = -ÊE The general displacement formula of the shuttle in this case can be expressed: 
Kk? 2 w K 
_ nee (L 3 2 = B A —Bt 
UETA (=) V E (cos(wt) +. sin(wt)) e | (14) 


considering the driving voltage with the sine wave function (Figure 4). 


Figure 4. V = Vo(1 + sin(Qt)) 


By setting: Fy = mse Vē, the moving force is given: 
0 
3 i 1 
Fas = Fo È + 2 sin( 2t) — + cos(29t)] (15) 


in this case, the deviation of dynamic model (4) is represented: 


3Fo 2Fo 


Vis = SE K-M2224C202 [(K = MQ?) sin( Nt) = CANcos(Qt)] 
et Fo = 2 3 
RETA [CK — 4M4 )cos(29t) + 2CQ sin( 2t)] 
+eF*IB, cos( wt) + B; sin(wt)] (16) 


according to initial condition: y,,(0) = Oandy,,(0) = 0, we have: 


B =F 2c K-4M02 3 

17 "0 {cman + 2(K-4M22)2+8C202 =} 

B, = “of 2cn? _ _2A(K-M2?) }+ £ B, 
w ((K—-4MR2)2+4C?02? (K-MR2)2+C22 w 


for the mathematical model (4) established aboveMy + Cy + Ky = F,(t). The control objective is to achieve 
the input signalF,(t) such that y(t) tracks y”°f (t). Due to the difficulties in finding the precise parameters 
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M,C,K ; it is assumed that the parameters M,C, Kare unknown and bounded by known constants Mo, Co, Ko, 
respectively. Let’s define that e(t) = {y"°f (t) — y(t)}, according to (4): 


d2 d d2 ref d ref 
MES +C% + Ke = (M +C + ky" — FO) 
that: 
2 
M3 +C&4+Ke = u(t) (17) 


2.,,7ref ref 
therefore, the proposed controller is chosen as F; (t) = (m 2 re +C ae + Ky" ) — u(t), where: 


de 
u(t) =D% (18) 
with D to be chosen as: 
2 
From (17) and (18), it implies that: M vo +(C-D) z +Ke=0 (19) 


hence, we only choose the parameter D such that the corresponding polynomial Mx? + (C — D)x + K = 0 
have all solutions being on the left side of complex plane. 


4. THE SIMULATION AND EVALUATION OF RESEARCH RESULTS 
The following simulations are performed to validate the tracking problem of output with the reference 
being the sinusoidal signal. The parameters of the proposed controller (18) are given as M = 1,C = 10, 


D=9,K =1. Two cases of sinusoidal reference signals are discussed as yı ref(t) = 2 sin (t+), 


Yoret(t) = 2 sin (t + z), The simulations will focus on the command tracking effect of the proposed controller 


(18). It is clearly seen that, although the references are depended on time but the tracking effectiveness is 
guaranteed (Figures 5 and 6). Moreover, the performance of tracking problem can be adjusted by choosing the 
parameters in (8) based on the roots of corresponding polynomial Mx? + (C — D)x +K = 0. 
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Figure 5. The result of the output with the reference being y,-¢(t) = 2 sin (t + z) 
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Response of Output 


Time (s) 


Figure 6. The result of the output with the reference being y,er(t) = 2 sin(t) 


5. CONCLUSION 

The article studies the dynamic model and control design of ELCA systems. Through the method of 
solving differential equation, the control scheme is implemented for ECLA in the presence of time varying 
reference. The tracking problem is analyzed through considering the solution of differential equation. By 
carrying out two simulation cases, it is concluded that the tracking problem is absolutely guaranteed with 
convergence rate to be depended on controller parameters. For future work, the problem of practical ECLA 
systems can be developed. 
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